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Abstract. The sky region including the Chandra position of SGR 1806-20 was monitored in the IR band during 2004, following 
its increased high energy bursting activity Observations were performed using NAOS-CONICA, the adaptive optics IR camera 
mounted on Yepun VLT, which provided images of unprecedented quality (FWHM better than O/l). After the 2004 December 
27th giant flare, the source position has been nailed by VLA observations of its radio counterpart, reducing the positional 
uncertainty to 0704. Using IR data from our monitoring campaign, we discovered the likely IR counterpart to SGR 1806-20 
\ based on positional coincidence with the Chandra and VLA uncertainty regions and flux variability of a factor of about 2 

correlated with that at higher energies. We compare our findings with other isolated neutron star classes thought to be related, 
' at some level, with SGRs. 
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^ ; 1. Introduction neutron stars (B ~ 10 14 -10 15 G; Duncan & Thompson 1992; 

^* ■ Thompson & Duncan 1995). The "magnetar" model is founded 

Soft Gamma-ray Repeaters (SGRs) were discovered in the sev- on two obse rvational facts: firstly, the rotational energy loss 



X 



enties through the detection of short (<ls), recurrent, and in- inferred from the SGR and Axp spin _ down is insufficient to 

^ ■ tense bursts of high energy emission peaked in the soft y rays. p0 wer their persistent X-ray luminosity of ~10 34 -10 36 ergs s" 1 ; 

Only four confirmed SGRs are known, three in the Galaxy and secondl y ; there is no evidence for a companion stars which 

one in the Large Magellanic Cloud (for a review see e.g. Woods could provide tne mass t0 power tne x _ ray emission through 

& Thompson 2004). The detection of a ~ 8 s periodicity in the accretion 
decaying tail of a very intense (~10 44 ergs) and long (several 

minutes) event, known as giant flare, from SGR 0526-66 on Bursting activity from SGR 1806-20 resumed at the end 

1979 March 5th (Mazets et al. 1979) suggested the associa- of 2003 displaying an increase in both the y-ray burst rate 

tion of SGRs with neutron stars. A small sample of peculiar and the hard X ~™y Persistent emission (Mereghetti et al. 

X-ray pulsars, namely the Anomalous X-ray Pulsars (AXPs) 2005a ) throughout 2004, and culminating with the giant flare 

has been proposed to be closely related to SGRs based on simi- of 27 } h December 2004 (Borkowski et al. 2004), during which 

lar properties, namely their period P (in the 5-12 s range), their ~ 10 er S s were released (for a distance of about lOkpc; 

period derivative P (lO^-lO" 13 s s" 1 range), and X-ray bursts Cameron et al. 2005; McClure-Griffiths & Gaensler 2005). Few 

(Kouveliotou et al. 1998; Kaspi et al. 2003; Gavriil et al. 2002). da y s after this event, SGR 1806-20 was observed and detected 

Both SGRs and AXPs have been proposed to be powered 
by the decay of strong magnetic fields that characterise these 

In this work we report on the results of an extended 

Send offprint requests to: gianluca@mporzio.astro.it Target of Opportunity (ToO) observational campaign on 

* The results reported in this Letter are based on observations car- SGR 1806-20 carried out during 2004 with the ESO VLT. In 

ried out at ESO, VLT, Chile (programs 072.D-0297, 073.D-0381 and particular, we report on the likely discovery of the IR counter- 

274.D-5018) part to SGR 1806-20 based on positional coincidence with the 



in the radio band for the first time, providing very accurate po- 
sitions (VLA; Cameron et al. 2005; Gaensler et al. 2005). 
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Table 1. Journal of the VLT NACO IR 2004 observations for 
the field of SGR 1806-20. 



2004 Start UT 


Filter Exposure FHWM Ks mag 


(MM DD / HH 


: MM) 




(s) 


(") 


for object A 


03 10 / 09 : 


00 


Ks 


1200 


0.14 


contaminated 


03 17 / 08 : 


14 


J 


1980 


0.21 




03 17 / 08 : 


53 


Ks 


1800 


0.10 


20.01 +0.14 


03 17 / 09 : 


33 


H 


1800 


0.14 




03 18 / 07 : 


35 


Ks 


3600 


0.09 


20.09 + 0.15 


06 13 / 03 : 


38 


Ks 


3600 


0.11 


20.26 + 0.26 


06 19 / 04 : 


51 


Ks 


6120 


0.11 


19.94 + 0.13 


08 10 / 02 : 


50 


Ks 


2640 


0.12 


contaminated 


08 11 / 03 : 


04 


Ks 


5160 


0.10 


19.48 + 0.12 


09 07 / 23 : 


56 


Ks 


5240 


0.09 


19.70 + 0.08 


10 05 / 23 : 


45 


Ks 


5160 


0.11 


19.32 + 0.16 



radio and Chandra positions and flux variability. (Preliminary 
results were reported in Israel et al. 2004, 2005a, 2005b, before 
and independently from Kosugi et al. 2005.) We briefly com- 
pare the IR emission properties of SGR 1806-20 with those of 
related objects. 

2. NAOS-CONICA observations at VLT 

The observations presented here were performed as part of an 
ESO Target of Opportunity program extending over October 
2003 - September 2004, and a Director Discretionary Time 
observation on October 2004. Observations were triggered fol- 
lowing the detection of intense y-ray bursts or during epochs 
of increased burst rate (Gotz et al. 2004; Golenetskii et al. 
2004a,2004b). The data were acquired at Yepun VLT with the 
Nasmyth Adaptive Optics System and the High Resolution 
Near IR Camera providing a pixel size of 07027 (NAOS- 
CONICA; see Table[T]for details). For all observations we used 
an exposure time of 40 s and a number of frames per image of 
3 with a random offset of 7" among images in order to perform 
background subtraction of the variable IR sky. VLT NACO sci- 
ence images were reduced based on the standard tools pro- 
vided by the ESO - Eclipse package (Devillard 1997). As a 
result of the presence of the Ks=8.9 magnitude object LBV 
1806-20 (see Eikenberry et al. 2004) close to the edge of the 
NACO field of view (FOV), artificial ring-like ghost structures 
were clearly detected in the image at coordinates R. A.= 1 8 h 08 m 
4031; Dec .=-20° 24' 41721 (equinox 2000), 13" away from 
the LBV 1806-20 position. In order to reduce as much as pos- 
sible the effects of contamination due to nearby objects, relative 
aperture and Point Spread Function (PSF) photometry was ob- 
tained within narrow annuli (about 1-1.5 FWHM depending 
on the seeing conditions), while the background was evaluated 
close to the object under analysis. Absolute photometry was de- 
rived by analysis of the best seeing frames. Finally, we cross- 
checked our absolute magnitudes by means of archival ISAAC 
data of the same region and about 100 isolated stars taken from 



the 2MASS catalog and within the instrument FOV: the results 
were in agreement to within 0.05 Ks magnitudes. 

In order to register the Chandra and VLA coordinates of 
SGR 1806-20 on our IR images, we obtained the image as- 
trometry by using the positions of about 10 stars selected from 
the 2MASS catalogs and within the -30" x 30" NACO FOV 
of final images. The residual in the fit was of 0706 in each co- 
ordinate, converting to ~ 071 once the 2MASS absolute accu- 
racy was included 1 . Fig.[T] shows the ~ 175x175 Ks band re- 
gion around the Chandra and VLA positions (l<x confidence 
level radius of Of! 3, Of! 04, respectively; Kaplan et al. 2002; 
Gaensler et al. 2005). However, given that the Gaensler et al. 
(2005) radio position refers to about 20 days after the giant 
flare of SGR 1806-20, and that the source from which is origi- 
nating the radio emission is moving at about 4 mas/day (Taylor 
et al. 2005), we also plot the VLA position obtained after 7 
days by Cameron et al. (2005; lcr radius if 071), corrected for 
about 30 mas in right ascension (following Taylor et al. 2005); 
this corresponds to a final lcr confidence level radius of Off 14. 
Source A, a relatively faint (Ks~20) object, at the sky position 
R.A = 18 h 08 m 39*337, Dec = -20° 24' 39785 (equinox 2000, 
90% uncertainty of 0706 ), is found to be consistent with the 
Chandra and VLA positional uncertainty circles superimposed 
on our IR astrometry-corrected frame. Objects B and C (-0723 
and Off 21 away from A, respectively; by looking at the contour 
lines, we note that object B might be the blend of two unre- 
solved objects) are only marginally consistent with the X-ray 
and radio positions, even though statistically plausible. Object 
A was not detected in the J and H images; 3<x upper limits of 
magnitude 21.2 and 19.5 were derived, respectively. We note 
that the SGR 1806-20 IR counterpart A (Ks magnitudes are 
listed in Tabled) plus objects B and C (Ks=19.07±0.04 and 
18.77+0.04, respectively) are all within a radius of -0725 from 
the VLA positions, and consistent with being unresolved com- 
ponents of candidate B in the IR images of Eikenberry et al. 
(2001;K=18.6±1.0). 

Light curves of the A, B and C objects marked in Fig. [T] are 
shown in Fig. [T] (right plot). Candidate A is the only one show- 
ing a clear brightening (a factor of ~2) in the IR flux between 
June and October 2004. Objects B and C show a fairly constant 
flux 2 . The upper panel of Fig.[T] (right plot) shows the closest 
reference star (176 away form the target) used for relative pho- 
tometry across Ks images: the object is constant to within the 
photometric uncertainties. We thus conclude that object A is 
variable. We checked for a similar variability also in the X-ray 
flux of SGR 1806-20. Both the XMM-Newton (Mereghetti et 
al. 2005b) and INTEGRAL (Mereghetti et al. 2005a) persis- 
tent fluxes of SGR 1806-20 showed an increase across the two 
semesters of 2004 by a factor of 1.94+°-™ and 1.7+^ in the 2- 
lOkeV and 20-100 keV bands, respectively 3 . During the same 



1 http://www.ipac.caltech.edu/2mass/releases/allsky/doc 

2 For the faintest component of blended object B we can reasonably 
exclude any IR variability, similar to that shown by A, for any Ks 
magnitude brighter than about 22.5 . 

3 for the 2-10 keV 2004 first semester flux we assumed that of 
October 2003, based on the unvaried INTEGRAL flux between 
October 2003 and February-April 2004. 
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Fig. 1. IR NACO Ks band image (left panel; 18th June 2004 dataset) close-up of the 175 x 175 portion of the sky around the l<x 
Chandra and VLA uncertainty circles (radius of 073, 07 '04 and 071, respectively) with the proposed counterpart marked with A. 
Coordinates are RA (hms) and Dec. (° ' "; equinox J2000). Isophotal contour lines are also drawn for clarity. Ks light curves 
(right panel) for a number of selected sources: the proposed counterpart to SGR 1806-20 (A) with two nearby objects (lower 
panel), and the closest reference (isolated and bright) star (upper panel) we used for relative photometry among Ks frames. Flux 
variability is clearly detected only for object A. 



time interval the NACO Ks flux increased by a factor of 2.4^, 
consistent with high energy flux variations. This further sup- 
ports the identification of object A as the correct IR counterpart 
of SGR 1806-20. 

Recently, independent from our work, the object A has been 
proposed as the IR counterpart to SGR 1806-20 (Kosugi et al. 
2005; their object B3). A comparison of their photometry with 
our we shows that nearly all the Ks magnitudes have an offset 
of about 0.2, with the important exception of objects A and C 
which are 1.6 and 0.6 Ks magnitudes brighter than the corre- 
sponding objects B3 and Bl in Kosugi et al. (2005), respec- 
tively. Even though we do not have a clear explanation for the 
observed differences, we note that in our images we did not see 
any evidence for (i) a brightening of objects B and C, and (ii) 
an increase of the local background around object A (based on 
our best datasets with FWHM <071), in contrast to Kosugi et 
al. (2005). An unusually high background level (regardless of 
its origin) may of course result in a flux underestimation of a 
source that lies in the same area. 



3. Discussion 

The deep and high spatial resolution NACO images allowed 
us to identify the likely IR counterpart of SGR 1806-20, and 
monitor its IR flux for seven months in 2004, during which 
an increase by a factor of ~2 was detected, correlated with 
the flux in the high energy bands. In fact, the IR flux of 
SGR 1806-20 was fairly constant until mid- June 2004, while 
it grew rapidly between June (Ks=20.01±0.14) and October 
2004 (Ks= 19.32+0. 16; l<x uncertainties; these values override 
the preliminary ones reported in Israel et al. 2005b). 



IR variability has been detected in nearly all AXPs with 
known IR counterpart. In particular, for IE 1048.1-5937, XTE 
J1810-197 and IE 2259+586, IR variability has been found, or 
suspected, to be correlated with the persistent X-ray emission 
(Israel et al. 2002; Rea et al. 2004; Tarn et al. 2004). Based on 
the NACO results we can conclude that the IR/X-ray correla- 
tion observed in AXPs also holds for SGR 1806-20. The total 
fluence of the IR enhancement between June and October 2004 
is about 10 41 ergs (we assumed Ay=29+2; see Eikenberry et al. 
2004), a factor of about 100 smaller than that in the 2-10 keV 
band. 

Based on the above reported findings we note that the 
SGR 1806-20 emission varies in a similar fashion (in terms of 
timescale and amplitude of variation) over more than five or- 
ders of magnitude in photon energy. The similar flux variation 
in the IR and X-ray bands suggests that the emission in the 
two bands has a similar, if not the same, origin. Moreover, it 
has become evident that X-ray flux enhancement of the per- 
sistent emission of SGRs is correlated with their burst rate, 
making it difficult to compare the fluxes among different SGRs 
without knowing their burst history (see Woods & Thompson 
2004). Tarn et al. (2004) argued that IR thermal surface emis- 
sion (within the magnetar model) is ruled out during the corre- 
lated X-ray/IR flux decay phases of IE 2259+586 (implausibly 
high implied brightness temperature), suggesting the magne- 
tospheric origin for the IR enhancement. Alternatively, the IR 
flux can be due to reradiation by material in the vicinity of the 
the pulsar. This model naturally predicts a correlation between 
the the IR and the X-ray flux (Perna, Hernquist & Narayan 
2000; Rea et al. 2004). 

This is the first time that the broad band energy properties 
of an SGR can be compared, over a similar energy band, with 
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Fig. 2. Broad band energy spectrum of SGR 1806-20 (tri- 
angles) and, as a comparison, the AXP IE 1841-045 (cir- 
cles) and the radio pulsar Vela (squares). In the case of 
SGR 1806-20 and IE 1841-045 high energy data are taken 
from XMM-Newton (Mereghetti et al. 2005b; referring to 6th 
October 2004 for SGR 1806-20), Chandra (for IE 1841-045) 
and INTEGRAL (Mereghetti et al. 2005a; 21st September - 
14th October 2004 for SGR 1806-20). Absorbed and unab- 
sorbed IR fluxes (A v =29±2, 5th October 2004 NACO obser- 
vation) are shown in the case of SGR 1806-20, unabsorbed 
(Av=13±l) IR fluxes are instead reported for the likely candi- 
date of IE 1841-045 (circles; Israel et al. 2005, in preparation). 
All the data for Vela are taken form literature (Kaspi, Roberts 
& Hardings 2004 and references therein). Solid curves (contin- 
uous, stepped and dot-stepped) are the unabsorbed fluxes for 
the black body plus power law model used to fit the high en- 
ergy part of spectra. 

those of other classes of isolated neutron stars, such as AXPs 
and radio pulsars. In Fig.[2]we show the "nearly simultaneous" 
broad band energy spectrum of SGR 1806-20 from the IR to y 
rays (high energy data are taken form Mereghetti et al. 2005a; 
see caption for details). The high energy part of the spectrum is 
clearly consistent with being non-thermal emission (a power- 
law model is generally used) from the source. We also plot the 
spectrum from the AXP IE 1841-045, for which 20-200 keV 
band data are available (Kuiper et al. 2004); a similar non- 
thermal component is displayed by the source. Non-thermal 
components are also seen in radio pulsars and modelled with 
power-law components (see Kaspi, Roberts & Harding 2004 
for a recent review). In some cases there is a smooth connec- 
tion between optical, X-rays and y-ray emission (Crab), while 
in other cases the extrapolation is plausible (Vela; see Fig.|2|). 
It is worth noting the similar flux ratios in the IR and hard X- 
ray bands for the three objects, and the significant difference 



of the characteristic temperature of thermal soft X-ray compo- 
nents between radio pulsars (<0.1 keV) and SGRs/AXPs (0.4- 
0.8 keV for a BB fit and 0.2-0.5 keV for a magnetic atmosphere 
fit, Perna et al. 2001), suggesting a significantly larger energy 
injection on the neutron star surface in "magnetar" candidates 
than in radio pulsars. 

Future detailed multi-wavelength observations campaigns 
of AXPs and SGRs will likely help clarifying the link between 
IR and high energy bands. Furthermore, the detection of the 
quiescent IR flux level of SGR 1806-20 will allow to compare 
the net energy released by the source in the IR and X-ray/y-ray 
bands during its bursting active phase. 
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